Cluster state generation using van der Waals and dipole-dipole interactions in optical 

lattices 
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We present a scalable method for generation of a cluster state for measurement-based quan- 
tum computing using van der Waals or dipole-dipole interactions between neutral atoms or polar 
molecules in an optical lattice. Nearest neighbor entanglement is accomplished by performing a 
phase gate using interaction of atoms in Rydberg states or molecules in large dipole moment states. 
All nearest neighbors are sequentially entangled in a finite number of operations, independent of 
the number of qubits, producing a ID cluster state. A universal 2D cluster state can be generated 
in several ms in a two-dimensional optical lattice by producing a series of ID cluster states in one 
lattice direction, followed by application of the entangling operations in another lattice direction. 
We discuss the viability of the scheme with Rb Rydberg atoms. 

PACS numbers: 



I. INTRODUCTION 

Entanglement plays a major role in quantum com- 
puting jl|, quantum communication [2[ and quantum 
metrology 0] . A special type of a multipartite entangled 
state, cluster state, represents a universal resource for 
the measurement-based quantum computing (MBQC) 
paradigm In measurement-based quantum comput- 
ing, quantum computations are carried out on the cluster 
state through individual qubit measurements in adap- 
tive bases (|0) ± e^> |1>)/V2. The primary advantage of 
MBQC over other register-based architectures rests in 
the fact that all interactions, required e.g. for two-qubit 
gates, are performed in the initialization stages of the re- 
source state. The computation may then be performed 
through simultaneous measurement of many individual 
qubits as warranted by the specific program being imple- 
mented. Any one- and two-qubit gate can be realized by 
appropriate measurements, making MBQC equivalent to 
the standard quantum circuit model [5(. Moreover, since 
the actual computation is done by local measurements, 
it can be faster compared to the equivalent gates in the 
circuit model. 

The cluster state is realized by preparing individual 
qubits in an eigenstate |+) = (|0) + |l))/\/2 of the 
Pauli spin operator a x , and letting the nearest neigh- 
bor qubits interact via (up to local operations) an Ising 

Hamiltonian H int — <?(*) X}( a ,a') * 2' ^~^t — during a 
time t so that JJ" g(t)dt = it. It clearly shows that 
the cluster state can be generated by applying a phase 
gate Upg = diag(l, 1,1,-1) between all nearest neighbor 
qubits. 

A cluster state was experimentally realized with pho- 
tons using linear optics, but this approach is difficult 
to scale to a large number of qubits. Schemes to gen- 
erate a cluster state using atom-cavity entanglement in 
the framework of cavity QED [7 L cold ions via phonon- 
mediated spin-spin interactions 8|, and distributed net- 



works of collectively excited atomic ensembles Q have 
also been proposed. A naturally highly scalable sys- 
tem is neutral atoms in an optical lattice, where a ID 
cluster state has been produced via ultracold s-wave col- 
lisions of atoms in a spin-dependent lattice [n|. Neu- 
tral atoms can also be entangled via long-range van der 
Waals (vdW) or dipole-dipole interaction when excited 
to Rydberg states. Atoms in Rydberg states can have 
huge dipole moments of several kiloDebyes and interact 
strongly with each other [ll[, which was suggested as a 
tool to produce two-qub it g ates in neutral-atom based 
quantum c omp uting: [12l 1 1 31 ] . Recently, in a series of ex- 
periments \lA\ Rydberg (dipole) blockade between two 
atoms has been observed, followed by the demonstrations 
of a blockade-assisted CNOT gate [15| and entanglement 
between two atoms [lj|. It is interesting to explore the 
possibility to use the strong interaction in Rydberg states 
to generate a cluster state. Moreover, it would offer a way 
to produce a cluster state with polar molecules, coupled 
via dipole-dipole interaction, for which collision-based in- 
teractions can result in undesirable inelastic or chemical 
reaction losses. For simplicity we primarily discuss in 
this work neutral atoms interacting in Rydberg states, 
the same scheme can be applied to molecules. 

We propose to realize a cluster state by applying a 
phase gate to all pairs of nearest neighbor atoms in an 
optical lattice. The phase gate can be realized using ei- 
ther direct or blockaded interaction in Rydberg states. 
This can be accomplished in a scalable way, i.e. in a fi- 
nite number of operations, not depending on the number 
of atoms in the lattice. In fact, to make pair excitation 
controllable and to minimize errors due to interactions 
of multiple atoms in Rydberg states, it already suffices 
to excite every other pair: (a) To produce a ID cluster 
state, four iterations are required and would be imple- 
mented as follows: using a periodic entangler, such as a 
standing wave, entangle positions 1 with 2, An + 1 with 
An + 2, etc.; in the second step, entangle positions 3 with 
4, An — 1 with An, etc.; in the third step entangle posi- 
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FIG. 1: (Color online) (a) An optical lattice described by 
the potential V = Vb cos 2 (fcx), along with an excitation 
pulse in the form of a standing wave with Rabi frequency 
Ql = Qlo cos 2 (kx/4 + <f>), where 4> — 7r /4; (b) Atoms are con- 
ditionally transferred from e.g. a qubit state |1) to a Rydberg 
state |r) by optical one or two-photon 7r-pulses. 

tions 2 and 3, An — 2 and 4n — 1, etc; in the last step 
entangle positions 4 and 5, An and An + 1 etc. (b) For 
2D cluster state generation, first the rows in x-direction 
are entangled, then the entangling operations are applied 
to columns in y-direction. This can be done using opti- 
cal standing wave excitation as is illustrated in Fig.[T^ 
in the case of a ID optical lattice. Pairs of atoms at the 
maxima of the standing wave intensity will get entangled, 
while those at the minima will not. Changing the phase 
of the standing wave fields, i.e. shifting the position of 
the intensity maxima and minima, will allow to perform 
the phase gate between all nearest neighbors, producing 
a cluster state. 

The paper is organized as follows. In Section II we 
describe how a ID cluster state can be realized using 
interaction in Rydberg states. In Section III we show 
how to generalize the technique used for the ID cluster 
state generation to produce a universal 2D cluster state. 
Finally, we discuss the main features of the scheme and 
conclude in Section IV. 



II. ID CLUSTER STATE GENERATION 

We start by analyzing how a ID cluster state can be 
generated by applying the phase gate between all neigh- 
boring atoms in a ID optical lattice. 

A. Phase gate without individual addressing of 
atoms in a pair (no dipole blockade) 

1. Gate analysis 

We assume that initially atoms are loaded into the 
ground motional state of an optical lattice described by 
a potential V(x) = Vq cos 2 (kx), shown in Fig.QJ,, which 
can be done using a superfluid-Mott insulator transi- 
tion. Two ground state hyperfine sublevels \F, tuf) en- 



code qubit states |0), |1). As a final preliminary step 
atoms are transferred into the |+) = (|0) + |l))/\/2 su- 
perposition by applying a n/2 pulse, resonant with the 
qubit transition. A phase gate between atoms in neigh- 
boring sites can be implemented using strong vdW or 
dipole-dipole interaction in Rydberg states. 

There are two ways to perform the phase gate depend- 
ing on whether atoms are individually addressable or not 
[12| . Standard Rydberg blockade requires individual ad- 
dressing of the atoms. Individual addressing in an optical 
lattice has been demonstrated recently using a tightly fo- 
cused laser beam and a spin-dependent lattice [l7| ■ How- 
ever, sequential application of the phase gate to each pair 
of atoms is not scalable to a large number of qubits. In 
a good scalable approach the cluster state has to be gen- 
erated in a finite number of operations which does not 
depend on the size of the system. Below we analyze the 
possibility to generate the ID cluster state assuming that 
atoms in each pair are not addressed individually. 

If one uses a standing wave excitation pulse with the 
Rabi frequency Hl — Qlo cos 2 (kx/A + 0), one has the 
maxima at e.g. even and minima at odd pairs sites (see 
Fig-Ht 1 )- The atoms at the maxima are transferred from 
the qubit state |1) to the Rydberg state \r) either directly 
or by a two-photon excitation, as shown in Fig.[Tp, be- 
low we assume two-photon excitation. The phase gate is 
then realized in the limit fi^ 3> Vint (Vint is the inter- 
action strength in the Rydberg state \rr)) as follows: 1) 
a 7r-pulse resonant to the |1) — |r) transition is applied 
simultaneously to both atoms, exciting each to the |r) 
state, since the shift of the doubly excited state \rr) is 
smaller than the Rabi frequency; 2) atoms in the \rr) 
state interact during time T- lnt and accumulate a ir phase 
shift Vi n tTi n t = 7r; 3) a second 7r-pulse deexcites atoms 
back to their original qubit states. The standing wave 
therefore allows to control the Rydberg excitation pat- 
tern, which would be more difficult if a spatially homo- 
geneous excitation pulse is used. We also note that vdW 
and dipole-dipole interactions of Rydberg atoms are long- 
range. As a result, when atoms are excited to |r), the 
interaction strength between atoms in neighboring pairs 
is comparable to the interaction strength between atoms 
within a pair, which will result in a phase error. The 
error is considerably reduced by exciting atoms only in 
every other pair using the standing wave. The interac- 
tion strength between atoms in closest excited even pairs 
is then smaller by a factor of <5VJnt/^int = 1/3 6 ~ 1CP 3 
for vdW and SV int /V int = 1/3 3 « 3.7 • 1CT 2 for dipole- 
dipole interaction than the interaction strength within a 
pair. Here Vint is the interaction strength between atoms 
in a pair and SVi n t is the interaction strength between 
atoms in closest exited pairs. 

Excitation to Rydberg states of alkali atoms is typ- 
ically a two-photon process via intermediate P1/2, P3/2 
states. Using the level scheme shown in Fig.[2ji we can 
write the Schrodinger equation for the amplitudes of the 
qubit states. If the qubit is initially in the |11) state, the 
corresponding system of equations is 



3 



.dan 2Vl\ 72^1^2 

i — ; — = : — «n : 

dt A A 



+ 



da 
i — 
dt 

I — — 

dt 



n 2 + n 2 



-an 



... 2f2 2 V2r!i^ 2 
(Vint - -£-)a r r ^ a H 



(1) 



where an, a+ and a rr are the amplitudes of the states 
1 1 1) , |+) = (|lr) + |rl))/v / 2 and \rr) (other states are far- 
detuned and not populated provided that A 3> Oi , ^2 ; 
A is the detuning from the intermediate state, £2^2 are 
the Rabi frequencies of the excitation pulses) . Assuming 
S7i = f?2 = £i and M n t <C £l 2 /A the solution of this 
system is given by 



|*) = e 2ia2 */ A 
|ll)-V2|+) + |< 



UD- \rr) 



-2iSl 2 t/A 



+ 



|ll)+V2|+) + |rr) c2 ^ t/A 



(2) 



which shows that for a pulse duration T such that 
fl 2 T/A = tt/2 the state evolves into |*> = |rr). Next, 
the excitation pulses are switched off and the atoms in- 
teract in the doubly excited state for time Tj nt so that 
Vint^lnt = 7T, and the state flips sign \*S?) = — \rr). Fi- 
nally, we apply the same 7r-pulse for time T, bringing the 
system into a state \ty) = — |11). 

If the initial state is 1 01) the system evolution is gov- 
erned by equations (similar for the 1 10) state) 



To proceed with cluster state generation the phase gate 
has now to be performed with odd pairs of atoms, i.e. 
atoms 3 and 4, An — 1 with An, etc. For that the phase 
4> of the excitation pulse is shifted by tt/2 to become 
(j) = 37T / 'A. At this point all atoms in even and odd pairs 
become entangled. As a next step the phase gate has to 
be applied to neighboring atoms, where one atom belongs 
to an even and another to an odd pair, i.e. to atoms 2 
and 3, An — 2 and An — 1, etc. and atoms 4 and 5, An 
and An + 1, etc. which can be done by setting <fi = tt/2 
and, finally, to <p = tt. This will result in the phase gate 
applied to all nearest neighbors. 

In Fig. [2b we show the level scheme of 87 Rb which 
we use as an example. A qubit can be encoded into 
|0) = |F= l,m F = 0) and \l) = \F = 2,m F = 0) states, 
providing long coherence lifetimes due to small sensitiv- 
ity to magnetic field fluctuations [l8|. A Rydberg state 
ns i/2 with \F = 2,m F = 2) can be used as \r), in this 
case atoms will interact via isotropic vdW Vint = Ce/r 6 
interaction. Atoms can be excited to \r) using two a + 
polarized pulses via intermediate pi/2 \F — 1,2, m F = 1) 
states. 



2. Phase gate errors 

The fidelity of the phase gate averaged over all initial 
two-qubit states is calculated in Appendix A in detail 
and is given by 



|(00|?7p G |00)| 2 - 
-(10|?7p G |10)| 2 



(01|t/ PG |01)| 2 + 
(11|L>p G |11)| 2 1, (6) 



= 1 v a aTa„7 ) ao1 a - ° 0r ' 

.da 0r /Ol tt\ \ n ± Q 2 fo , 

= U aTa^ ) aor - ~ ° 01 ' (3) 

which in the case A 3> Ahf gives the solution 
|^ = e 2 l n 2 t/A / [01) - \0r) c _ iQ 2 t/A | 

m) + \°r) j*t,^ t (4) 

where Ahf is the hyperfine splitting of the atomic ground 
state. One can see that for fl 2 T/A = tt the system re- 
turns to the state |vfr) = - |01). 

Finally, if the system is initially in the 1 00) state, the 
wavefunction evolves as 

|*) = e 2 ^/(A+A ht ) | 0Q) ; (5) 

bringing the system into |^) = 1 00) for Q 2 T/A = tt 
if A ^> Ahf. As a result, the phase gate \e1e2) — > 
_ e »7r(i-ei)(i-£ a ) | ei£2 ) is implemented. 



where the state after the imperfect phase gate Upg is 
compared to the state after the ideal gate |11) — > — |11), 
|01) - |01), 1 10) - 1 10), 1 00) ->■ |00). In Appendix 
A we calculate several types of intrinsic errors and find 
that the main errors of this type of phase gate are caused 
by the finite width of the ground motional state wave- 
function of each atom in a lattice site. This results in a 
finite spread of the Rabi frequency of the standing wave 
excitation pulse. Other significant errors are due to the 
finite ratio of the interaction strength to the two-photon 
Rabi frequency and the decay of population in Rydberg 
states. The corresponding gate fidelity is 

F=l-e=l- 2TT 2 ({Sfl/(n)) 2 ) 

-V 2 J(8(n 2 /A) 2 ) - 2TTJ/V int - ttj/^/A), 

where e is the corresponding gate error, and averaging 
is over the motional ground state wavefunction. In the 
fidelity calculation equal Rabi frequencies of the pulses 
£2i = O2 = f2 were assumed and a decay of a Rydberg 
state with the rate 7 was introduced. 

In order to find the error due to the finite width of 
the ground state motional wavefunction we approximate 
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the potential at each site as har monic with the oscil- 
lation frequency u> — ky /, 2Vo/m and the correspond- 
ing wavefunction width a = (Er/Vo) 1 ^ 4 /k « 0.316/fc 
for Vo = 100 E a, which we use for an estimate. Here 
Er = h 2 k 2 /2m is the atomic recoil energy, m is the 
atomic mass. We assume the Rabi frequencies of the 
excitation pulses Six = Sl 2 = SI ~ cos(fcir/4 + 7r/4), where 
x varies around potential minima of the n th even pair of 
sites x n — 57r/2(77r/2) + Ann with the Gaussian proba- 
bility distribution p(x) = exp(— (x — x n ) 2 / a 2 ) / ^Jita. 
The contribution of the first error term is 

eovar = 27r 2 (( ( 5r!/(r!>) 2 > = 

= 7T 2 (fca/4) 2 tan 2 (kx n /4 + ir/4) 1.05 • 1(T 2 . 

Choosing V int = 3 MHz and 2 /A = 30 MHz, the 
second, third and fourth error terms are ei mp exc = 

l/8(Mnt/^ 2 /A) 2 » 1.25- 10- 3 , CRydb decay 1 = « 

6.7- 10- 4 and e Rydb dccay2 = 7r 7 /(f! 2 /A) « 3-10- 5 for the 
Rydberg state lifetime 500 /its, i.e. the total error is of the 
order of e « 1.25 TO -2 . The high value of the two-photon 
Rabi frequency f2 2 /A = 30 MHz can be achieved using a 
two-photon excitation path — > 6pi/2 ~^ ns i/2 with 
422 and 1004 nm lasers, respectively, due to the larger 
dipole moment of the Qpi/2 ns\/2 transition [l9l ]. al- 
lowing to reach a high Rabi frequency Q2 ■ At the same 
time the decay time of the 6pi/2 state (125 ns) allows to 
minimize the error due to the decay of the intermediate 
state. The probability of the decay during the gate is 
Psc = TTTepirt/A- Choosing A = 40 GHz the probability 
of the intermediate state decay is p se ~ 10 -4 . 

There is also an error due to the undesirable interaction 
of atoms belonging to different excited pairs, e^if pairs = 
(3^ 2 /16) (1/8+ 19tt 2 /256) (SV int /V int )(V int /(n 2 /A)) 3 
(see Appendix A, section A). Assuming that inter- 
acting atoms belong to the closest excited pairs, 
we have SV int /V int = 1/3 6 w 10~ 3 for vdW and 
SV int /V int = 1/3 3 3.7 • 10" 2 for dipole-dipole inter- 
action. As a result, the error is £dif pairs ~ 1-6 • 10~ 6 
for vdW and edif pairs ~ 5.9 • 10 -5 for dipole-dipole 
interaction. 

This analysis assumes that the excitation from |1) 
to |r) and back is adiabatic with respect to the lat- 
tice motional frequency and there is no uncertainty in 
the interaction strength, analyzed in Appendix A, sec- 
tion C. If the excitation-deexcitation is non-adiabatic, 
there is an additional gate error due to the uncertainty 
in the interaction strength e non a diab = ( 7r2 /4)(a/-fi > ) 2 — 
(7r 2 /4)(a/i?) 4 (Vint/w) 2 , where R is the distance between 
nearest neighbors in the lattice. This error was derived 
assuming V5 nt < u>, which is not the case in our esti- 
mates. If Vj n t > uj the error will be higher, and we use 
the first term e non adiab = (tt 2 /4)(a/ R) 2 s» 2.5 • 10 -2 as a 
lower bound for this type of error. The error derivation 
assumes that the trapping potential for the ground and 
Rydberg state is the same, providing the same motional 
frequency. This can be achieved in a blue-detuned optical 
lattice at a " magic" wavelength, at which ground and Ry- 
dberg state polarizabilities are equal [20]. On the other 




(a) (b) 

FIG. 2: (Color online) (a) Two-photon excitation from the 
|11) to the \rr) state in alkalis; (b) Level scheme of 87 Rb 
showing qubit |0), |1) and Rydberg \r) states. 



hand, the error can be avoided if the trapping lattice is 
switched off when atoms are excited to Rydberg states 
and the gate pulses are much shorter than the motional 
period. In this case atoms move only a small fraction 
of the ground motional state width during the gate and 
will be recaptuted in the ground motional state once the 
trapping lattice is switched back on. 

Finally, there is a finite probability to excite atoms in 
"inactive" lattice sites, sitting in minima of the standing 
wave excitation pulse. The probability that the pair of 
atoms stays in the initial state after the gate, averaged 
over four possible initial states, is calculated in Appendix 
A, section E. For a two-photon excitation the averaged 
probability (P) w 0.75. 

The error analysis shows that first, vdW interaction re- 
sults in smaller errors caused by the interaction of atoms 
in different excited pairs compared to the dipole-dipole 
one. Second, the excitation to Rydberg states has to be 
adiabatic to avoid errors due to the uncertainty of the in- 
teraction strength. On the other hand, these errors could 
be avoided altogether if excitation to the \rr) state is not 
required. Next we therefore analyze the phase gate based 
on Rydberg (dipole) blockade, where only one atom is ex- 
cited to a Rydberg state, provided individual addressing 
of atoms in a pair is possible. 



B. Phase gate with "individual" addressing: dipole 
blockade 

1. Gate analysis 

If individual addressing of atoms in a pair is possible, 
one can use standard dipole blockade [l2j , and the phase 
gate is realized in the limit fl l -C Vi n t as follows: 1) 
a control atom is excited from one of the qubit states, 
e.g. a qubit state |1) to a Rydberg state \r) by a ir- 
pulse; 2) a 2-7r-pulse of the same frequency is applied to a 
target atom. Dipole-dipole or vdW interaction shifts the 
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energy of the doubly excited state \rr), as a result the 
27r-pulse has no effect on the target atom (the excitation 
is blockaded) since the Rabi frequency of the pulse is 
much smaller than the energy shift; 3) the control atom 
is deexcited back to its original qubit state. We stress 
that we do not require individual addressing of atoms 
in the lattice, we only need atoms in each pair to be 
separately addressable, this is why we call it " individual" 
addressing. 

Individual addressing in a pair can be realized using a 
polarization gradient optical lattice plj . where not only 
the intensity of the lattice field but also its polarization 
changes periodically in space. A lattice of double-wells 
was demonstrated in [22| with the polarization being lin- 
ear in one site and elliptical in the neighboring site of each 
well. The elliptical polarization results in a non-zero con- 
tribution from a vector part of polarizability, producing 
a state-dependent shift ~ ot v (E L x El)F, where a v is 
the vector part of the polarizability, F is the atomic an- 
gular momentum, and El is the positive frequency part 
of the total electric field E = £^e~ iwt + c.c. The contri- 
bution from the vector part can, therefore, be viewed as 
a fictitious magnetic field Bg ct ~ a v {E* L x El). 

Atoms can be loaded into a double- well optical lattice, 
formed by two standing wave fields of different polariza- 
tions. Let us consider the lattice field 



E L = E e y e lkx+l 

+iE 1 e z e 2ikx 



E Q e y e- tkx - 14 



iE\e z e 



-2ikx 



This lattice can be produced by two pairs of laser beams, 
intersecting at tt and 7r/3 angles (23|. The corresponding 
lattice potential is a sum of scalar and vector parts V = 
-a s \E L \ 2 /4:+ia v (E* L xEL)F/4 = V s +V v , where a s is the 
scalar polarizability. The scalar part V s = Vbcos 2 (A;a; + 
4>) + Vi cos 2 2kx (where Vq — —a s El and V\ = —a s E 2 ), 
represents a double- well lattice, with a spatial period a — 
n/k; the minima of the n th double- well are at (kx m i n ) n = 
(arcsin(K) /4Vi ) /2 + 7r/4+7rn, - arcsin(F /4Vi)/2 + tt/4+ 
7r(2n + l)/2). The vector part is given by 



V v — ~2a v \/VoVi cos(2kx)cos(kx) F x /a s , 



(7) 



and at the minima of a double-well V v = 
±a v y/V Vi(l - V /4Vi)/2V F x /2Via a . Here F* is 
the projection of the total angular momentum F on 
the direction of the fictitious magnetic field Bfi ct (x- 
direction), i.e. F x = mp] the ± signs refer to the left and 
right sites of the well. The scalar part of the potential 
is typically larger than the vector part, since the vector 
polarizability is about an order of magnitude smaller 
than the scalar part for alkalis. Assuming a v /a s ~ 0.1, 
Vo = V\ and Vq = lOOE^ the vector shift at the well's 
minima is V v « ±3ErF x . We can encode a qubit into the 
states |0) = \F= l,m F = -1), |1) = \F = 2,m F = 1), 
this qubit transition is insensitive to magnetic field 
fluctuations at a bias magnetic field B^ ms = 3.23 G 
(l8| . As a result, if a magnetic field is applied, the |0) 



and |1) states will shift by the same amount preserving 
the qubit frequency. The V v term, acting as a fictitious 
magnetic field, will then shift states |0), |1) in the right 
site with respect to the qubit states in the left site of 
the well [24j (see Fig-EJa), allowing to selectively excite 
only one atom in each pair, while having another atom 
unaffected by the excitation pulse. For 87 Rb Er = 3.5 
kHz at the lattice wavelength A = 810 nm, and the shifts 
for the |0), |1) states are V v ~ ±10 kHz. Larger shifts 



24] can be realized with 



of the order of 50 - 100 kHz 
optimized parameters of the lattice. 

If the two-photon excitation Rabi frequency is much 
smaller than the vector shift A voc = 2\V V \, the atoms in 
the left and right sites can be excited to the |r) state 
selectively. The phase gate then can be realized using 
dipole blockade as: 1) a 7r-pulse resonant to the |1) — |r) 
transition of a control atom is applied, exciting it to the 
|r) state; a target atom is off-resonant by A voc , which is 
the relative shift of the qubit state |1) in left and right 
site, and is not excited; 2) a 27r-pulse resonant to the 
|1) — \r) transition of the target atom is applied. The 
doubly excited \rr) state is shifted by the large interac- 
tion energy Vint and is not populated. The control atom 
is off-resonant by A voc and is not affected; 3) finally, a 
7r-pulse resonant to the control atom brings it back to its 
original qubit state. These steps are analyzed in detail 
in Appendix B and are summarized in the table below, 
which shows the evolution of the two-qubit states: 



|00) ^ 


-|00) 


27Tf 


- 100) ^ 


|oo> , 


|01) ^ 


-|01) 


2tt, 


|01) ^ 


-|oi), 


|10) ^ 


-i\r0) 




-ie %e \r0) ^ 


-e ie |10), 


111) ^ 


-i\rl) 


2ti-( 


-ie l8 \rl) ^ 


-e* e |ll), 



(8) 

where 6 = 7rA vcc /fi 2 / A, the first qubit corresponds to 
the control and the second to the target atom. One can 
see that the phase gate can be realized provided that 
9 = 27m. 



2. Phase gate errors 

The averaged fidelity of the phase gate with dipole 
blockade was calculated in Appendix B and is given by 

F=l-e=l- 2Tr 2 ({Sn/(n)) 2 ) - 77r 7 /(4fi 2 /A) 

-(fi 2 /A) 2 /(2A 2 0C ). 

The error terms are due to the spread of the two-photon 
Rabi frequency experienced by each atom 



envar = 27r 2 ((cftVr!) 2 H 
0.15, 



1 + V(T+ Vo/AV^/2 

decay of Rydberg states e Ryd b decay = 7TTj/(An 2 /A) rj 
4.38 • 10 -2 , and the imperfect frequency selectivity 
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FIG. 3: (Color online) (a) Double-well lattice described by 
the potential V = Vo cos 2 kx + V\ cos 2 2kx (Vi = Vb) along 
with a standing wave excitation pulse with Rabi frequency 
Q, L = fl L0 cos 2 (kx/2 + 7r/4); (b) Shift of qubit states |0) = 
\F = l,mj7 = -1} and |1) = \F = 2,m F = 1} in left and right 
sites of a double-well in a polarization gradient lattice de- 
scribed in the text. 



between the control and target atoms €i mp block = 
(17 2 /A) 2 /(2|A VCC | 2 ) « 0.02. Here we assumed Vint > 
A vcc » ft 2 / A, v n = V x , a = (E B jV )^ 4 /k « 0.316/* 
for Vo = IOOEr, I/7 = 500 ^s, the two-photon Rabi fre- 
quency ft 2 /A = 40 kHz, and the vector shift of the qubit 



states A v 



200 kHz. It shows that the total error of 



the phase gate in this case is e ~ 0.21 and the phase gate 
time is T PG = 2?r/ft 2 /A pa 25 /xs. 

Let us also discuss the error due to undesirable exci- 
tation of atoms in "inactive" wells, where the Rabi fre- 
quency is close to a minimum. It was calculated in detail 
in Appendix B, section D, and the probability to find the 
pair of atoms in the initial state after the gate, averaged 
over all four initial states, is (P) ss 0.998 in the case of 
two-photon excitation to Rydberg states. 

By changing the phase of the standing wave so that in- 
tensity maxima shift to odd double-wells the phase gate 
can be realized in odd pairs. At this stage all pairs in 
double-wells are entangled. In the next subsection we 
describe how the phase gate can be realized with neigh- 
boring atoms belonging to different double-wells. 



3. Lattice manipulation 

In the previous subsection we showed how a phase gate 
can be applied to pairs of atoms in double- wells. To 
proceed with the cluster-state generation the phase gate 
has to be realized with the neighboring atoms in different 
double- wells, i.e. between each atom in a right site of 
the n th double- well and an atom in a left site of the (n + 
l) th double-well. The phase gate operations described in 
the previous subsection can be applied if the atoms are 
brought to the same double-well. This can be achieved 
by adiabatically manipulating the lattice in the following 
way: (i) decreasing Vq (Fig. Hp, left panel) which raises 
the barrier in each double-well, (ii) ramping the phase <f> 
from to 7r/2 (Fig.0t, left panel), (iii) finally, increasing 




FIG. 4: (Color online) Left panel: (a) Double-well lattice 
described by the potential V = Vo cos 2 (kx + (j>) + Vi cos 2 (2kx), 
Vx = Vo, <t> = 0; (b) Vo is ramped down to Vb = 0.2Vx; (c) The 
lattice phase is shifted from to <f> = 7r/2, Vb = 0.2Vb; (d) Vb 
is ramped up to Vo = Vx', Right panel: populations of the four 
lowest motional states during the Vb and <j> manipulations. 



Vo back to its initial value. As a result, the atoms that 
were in the right and left sites of neighboring double- wells 
end up in the same double- well (Fig.HJi, left panel). 

The lattice manipulation has to be adiabatic to avoid 
motional excitation of atoms. Modeling the evolution of 
the atomic motional state during the lattice manipula- 
tion requires calculation of Bloch bands and Bloch func- 
tions for every configuration of the lattice. These are 
the eigenstates and eigenfunctions of the single-particle 
atomic Hamiltonian, which can be found by solving the 
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Schrodinger equation with the potential V(x,t). The 
potential V(x, t) is periodic, and can be written as a 
discrete Fourier sum, containing terms with wavevectors 
k x = ±fc,±2fc, in the following way: 



V(x,t) = (Vo + V^/2 + Voie 



2ikx+2i<p 



+ e 



)/2- 



+V 1 (e 



4/k.L 



+ e 



—4ikx 



)/2, 



allowing one to write the solutions of the Schrodinger 
equation in the form of Bloch functions ipq(x) — 
e lqx u q (x). Here q is the quasi- momentum (restricted to 
the first Brillouin zone), n is the band index, and 



UniX) 



=-iV n 



2inkx 



(9) 



with iV max a suitable cutoff number. The resulting sys- 
tem of equations for the c„ coefficients and eigenenergies 



h 2 {q + 2nkf 
2m 
Vi 



+y(c„ +2 +c„- 2 ) = (£ (n) (q) 



V Q + V X 



is solved numerically by truncating the sum in Eq.([9"|) 
at some iV max providing a necessary precision for the 
eigenenergies. We used N m ax = 10 to calculate the low- 
est Bloch energies and functions. Different configurations 
of the double-well potential and the corresponding band 
structures are shown in Fig. [5] The left panel demon- 
strates Bloch bands corresponding to the three stages 
of lattice manipulation, while in the right panel of the 
figure the energies of the q = eigenstates of several low- 
est bands are shown. The energies of the lowest bands 
weakly depend on q, as can be seen from Fig.[5fa) ) (c) ) (e), 
and the energies and Bloch wavefunctions corresponding 
to q = can be used. Given the Bloch functions tpq n \ 
which are delocalized over the entire lattice, one can con- 
struct Wannier functions which are localized at lattice 
sites Xi 



w 



(n) 



1 



~ lqx ^l n \x). 



(10) 




FIG. 5: (Color online) Left panel: Bloch bands of the double- 
well lattice V = Vo cos 2 (kx + (j>) + V 1 cos 1 (2kx), (a) Vi = Vo = 
IOOSh, <j> = 0; (c) Vo = 0.2Fi, = 0; (e) Vo = 0.2Vo, <f> = tt/2; 
Right panel: Double-well potential along with lowest Bloch 
energies, where Vo and 4> m (b), (d) and (f) are the same as 
in the left panel for (a), (c) and (e). 




In a double- well lattice the two lowest bands are sepa- 
rated in energy by much less (< O.IEr) than the typical 
motional excitation energy ~ y/AVoER. As a result, even 
at ultracold temperatures both bands are going to be 
populated. In this situation to obtain a correct descrip- 
tion of the system evolution generalized Wannier func- 
tions are introduced [25], which are superpositions of the 
Wannier functions of different energy bands. In our case 
the Bloch functions for the first two bands are symmetric 
(ground band) and anti-symmetric (second band) around 
the center of a double- well (shown in Fig.[S^,) . Combining 
the Wannier functions corresponding to the two bands as 
ipL,i = (ipi.i - ^2,i)/V2 and ip R i = (ipi ti + ip 2 ,i)/V^, one 



FIG. 6: (Color online) (a) Bloch functions tpi (green dot- 
dashed curve) and ip2 (red dashed curve) of the first and sec- 
ond Bloch bands along with the lattice potential 10 _2 V; (b) 
Wannier functions centered in the left = — 4'2,i)/V2 
and right tpRj = (ipi.i + ip2,i)/V2 sites of a double-well. 



can obtain generalized Wannier functions localized in the 
left and right well, respectively (see Fig.|6p). 

The lattice parameters have to be changed adiabat- 
ically to avoid undesirable motional excitations of the 
atoms, i.e., slow compared to the lattice motional ener- 
gies. We checked the adiabaticity of lattice manipulation 
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by numerical modeling of the system evolution. We took 
into account first four Bloch states to simplify the anal- 
ysis and assumed that initially atoms are in the ground 
Bloch state. As we already mentioned, the lowest Bloch 
bands weakly depend on the quasi-momentum q, and, as 
a result, one can approximate the corresponding Wannier 
functions by the Bloch functions ipL,i and ipR,ii restricted 
to a single site. The population evolution is shown in the 
right panel of Fig. [4] along with the time-dependence of 
the amplitude Vo and phase <$>. By optimizing the ma- 
nipulation time we found that 99.92% of the population 
stays in the ground Bloch band. We note that this value 
can be further increased by using more complex (opti- 
mized) functions Vo (t) and <j>{t) . The total time required 
for the manipulation is ~ 600 /is. 

We can now estimate the time required to generate 
a ID cluster state. In the case considered in Section 
II A, when atoms are not addressed individually, the 
phase gate has to be performed adiabatically with the 
gate time Tpc 3> 1 /is, as a result, the excitation- 

deexcitation pulses have to be longer than 10 /is. We can 
assume the phase gate duration Tpq « 20 ^s. The result- 
ing time of cluster state generation, including four phase 
gate sequences is then T\d sa 80 /is. If atoms in a pair 
are individually addressable, the gate duration was found 
to be w 25 /is in Section II B. The total time including 
the lattice manipulation is then T\d ps 700 /is. 



III. 2D CLUSTER STATE GENERATION 

The scheme described in the previous section can be 
extended to generate a 2D cluster state, which is required 
for universal quantum computation. First, atoms can 
be loaded in a 2D lattice with V = Vocos 2 (fcix + <f>) + 
Vi cos 2 (2fcix) + V% cos 2 (2fc 2 y), which produces a regular 
lattice with a period it /k 2 in y direction and a regular or 
double-well lattice in x direction. The y-lattice period is 
assumed sufficiently large so that atoms in neighboring 
x-chains do not interact when excited to Rydberg states. 
In this way, following the steps of Sections II A or II B, 
we can produce a series of ID cluster states in the x di- 
rection. As a next step we adiabatically reduce V\ to get 
a regular lattice in the x-direction with a period n/2ki 
(in the case of a double- well lattice), followed by stretch- 
ing of the x-lattice. Next, the y-lattice period is reduced 
to bring the x-chains closer and the regular or double- 
well lattice V = V% cos 2 (2fc 2 y + 4>) + V 3 cos 2 (fc 2 y + 4>) is 
switched on in the y direction. Next, the entanglement 
operations of Sections II A or II B can be repeated for ID 
chains in the y direction, producing a 2D cluster state. 

Let us estimate the time required to generate the 2D 
state. As we showed in the end of Section II C a ID clus- 
ter state in the x-direction requires Tid ~ 80 /is with- 
out and Tin ~ 700 /is with individual addressing in a 
pair to be realized. Adiabatic reduction of Vo will take 
sa 250 /is, as can be seen from Fig.|U[e) in the right panel. 
Next, we need to estimate the time required to stretch 



the x-lattice. The stretch has to be adiabatic for atoms 
to remain in the ground motional state of the lattice. 
We modeled the stretch of the lattice V = ViCOS 2 (fcx) 
so that the k vector was adiabatically changed from 2ki 
to 0.4fci. The period of the lattice 5n/2ki at the end 
of the stretch is comparable to the distance 37r/2fci be- 
tween atoms which belong to every other excited pair in 
the double- well. As a result, the error due to the inter- 
action of atoms in neighboring x-chains is smaller than 
the error due to the interaction of atoms in neighboring 
excited pairs in the same chain. 

We modeled the evolution of the system by calculating 
the Bloch bands and Bloch functions, and used the latter 
to construct single-site Wannier functions of the lattice 
V(x, t) during the stretch. We again took into account 
four lowest Bloch bands, assumed that initially the popu- 
lation was in the ground Bloch band, and monitored the 
excitation to higher-energy bands. We found that the 
stretch can be performed rather fast in ~ 16 /En ~ 730 
/is, where En = 3.5 kHz for 87 Rb, while retaining 99.55% 
of the population in the ground Bloch band. The popu- 
lations of the lowest four bands during the lattice stretch 
are shown in Fig.[7J 

Finally, we can estimate the total time required to gen- 
erate the 2D cluster state in the blockaded case. The time 
required for ID state generation T\r> «a 700 /is; ramping 
down the Vo lattice takes rj 250 /is and the lattice stretch 
in the x-direction takes sa 730 /is; the total time is sa 1.7 
ms. Next, the lattice period in the //-direction has to 
be adiabatically reduced, which will similarly take pa 730 
US, followed by ramping up the V3 lattice in pa 250 /is, 
and a sequence of phase gate operations applied to the 
//-chains, will take Tid ~ 700 /is. Therefore, the total 
time of the 2D cluster state generation in the blockaded 
case is T 2 d ~ 3.4 ms. The time required for 2D clus- 
ter state generation in the no-blockade case can be found 
summing the time required to produce ID states in x- 
direction, Tid pa 80 /is, the time required to stretch the 
lattice in x-direction and shrink it in y direction pa 2 x 730 
/ts=1.46 ms, followed by application of the phase gate se- 
quence in y-direction, requiring time T\d ps 80 /is. The 
total time needed to produce a 2D cluster state is there- 
fore T 2 d * 1-62 ms. These times are much smaller than 
the qubit coherence times ~ 1 s [l8j , achievable in optical 
lattices. 



IV. DISCUSSION AND CONCLUSIONS 

In Sections II. A. 2 and II. B. 2 we showed that the phase 
gate between two neighboring atoms can be realized with 
an error e pa 1.25 • 10 -2 and e «a 0.21 without and with 
dipole blockade in the pair, respectively, if two-photon 
excitation to Rydberg states is used. The major contribu- 
tion to the error in both cases comes from the variation of 
the two-photon Rabi frequency experienced by each atom 
due to the spatial variation of the excitation held. This 

error can be reduced to en var = (7r 2 /2) ( (5Q,/(£l)) 2 V i.e. 
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FIG. 7: (a)-(d) Populations of the four lowest Bloch bands 
during the lattice stretch; (e) the k vector is linearly reduced 
from 2fci to QAki . 



four times, if a single-photon excitation is used instead 
of a two-photon one. Then, en var ps 2.63 ■ 10~ 3 and 
en var ~ 3.75 • 10 -2 in the cases without and with dipole 
blockade, respectively. It will reduce the total gate error 
to e ps 4.3 • 10~ 3 without and e pa 0.16 with dipole block- 
ade. The error of the blockaded phase gate can be further 
reduced by increasing the relative qubit energy shift A ve c 
between left and right sites in a double well. Namely, the 
error term £i mp block = (l/2)(f2/A vec ) 2 (here one-photon 
excitation is assumed) can be reduced by increasing A vec . 
One possible way to increase A vec is to use hyperfine 
states \F,mp) with large \rriF | for qubit encoding, since 
A vcc ~ nip. For example, if Cs \F = &,mp=A) and 
\F — 3,mp = 3) states are used as qubit states |1) and 
|0), this will allow the increase up to A VO c ~ 200 kHz x 
(nip =4) ~ 800 kHz. The A voc can be increased even 
further using atoms with large hyperfine quantum num- 
bers F, such as the rare-earth Ho, having 4 < F < 11 
in the ground 4/ 11 6s 2 ( 4 / 1 5/ 2 ) state H3|, which will al- 
low the increase up to A vec ~ 1 MHz. Choosing one- 
photon Rabi frequency il ~ 100 kHz, the blockaded gate 
errors become e imp block = ^ 2 /(2A 2 0C ) ps 5 • 10~ 3 and 
£Rydb decay = 77T7/4S1 ps 1.75 • 10 -2 . The error due to 
Rabi frequency variation en var 



(7r 2 /4)(fca) 2 



-^/(l + V„/iV 1 )/2 



is harder to reduce, since 

l+y/(l+V /4,V 1 )/2 

it pretty much depends only on the ratio Vo/Vi, which 
cannot be increased much beyond Vo ~ V\ in order to 
have motional states localized in left and right sites of 
a double-well. Using a very deep optical lattice with 
Vo = 200-Er and Vq = 2Vi, which still results in mo- 
tional states localized in left and right wells, allows to 
reduce the error to en var ~ 1-13 • 10~ 2 . The transition 
between \F = A,mp = 4) and \F — 3,mp = 3) states of 
Cs is sensitive to magnetic field at all field values, mean- 
ing that the qubit will experience decoherence in these 
states. To reduce its effect one can keep the qubit in field 
insensitive \F = 4, mp — 0) and \F = 3, mp — 0) states 
during storage and lattice manipulation time and trans- 
fer it to the \F = 4, mp — 4), \F = 3, mp — 3) only when 
the phase gate is applied. This will result in an addi- 
tional phase gate error due to magnetic field fluctuations 
£mf fluct ~ (<5wTpg) 2 , where 5u ps g e ^sSB is the fluctu- 
ation of the qubit transition frequency due to the fluc- 
tuating magnetic field SB and g e is the electron's gy- 
romagnetic ratio. In a recent study [28[ suppression of 
magnetic field fluctuations down to SB ~ 50 /iG was 
demonstrated, resulting in Suj ~ 10 3 s _1 . One phase 
gate operation requires Tpo ~ 25 /xs in our setup, result- 
ing in the dephasing error per gate cmf fluct ~ 6 • 10 -4 . 
As a result, the total gate error in the blockaded case can 
be reduced to e ps 3.44 • 10~ 2 . 

At the same time there are errors due to the excita- 
tion of atoms in "inactive" latice sites, situated at the 
minima of the standing wave. With a two-photon ex- 
citation, the probability to find the pair of "inactive" 
atoms in the initial two-qubit state after the gate, av- 
eraged over all four initial states, is (P) ps 0.75 and 
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(P) ps 0.998 in the cases without and with blockade, re- 
spectively. The probability in the former case can be in- 
creased using multi-photon, e.g. four-photon excitation 
to the Rydberg states. With a four-photon exitation, 
the probability increases up to (P) s» 0.994 in the no- 
blockade case. At the same time the error due to the 
variation of the Rabi frequency increases four times up 

to e n var = 8tt 2 ((6n/(Q,)) 2 ^ 4.2 • 10~ 2 , giving the to- 
tal gate error e sa 4.4 • 10~ 2 in the no-blockade case. On 
the other hand, single-photon excitation which reduces 
the gate error in the blockaded case, leads to the reduced 
probability to leave atoms in "inactive" sites unaffected 
(P) « 0.87. 

We can conclude from this analysis that in no-blockade 
case the multi-photon excitation, e.g. four-photon, leads 
to the optimal combination of the phase gate error (e pa 
4.4- 10 -2 ) and the error due to the excitation of atoms in 
"inactive" lattice sites (ei nac t cxc = 1 — (P) ~ 6.5 • 10 -3 ). 
At the same time, in the blockaded case one-photon ex- 
citation seems to provide the optimal error combination: 
phase gate error e 3.44 • 10 -2 which can be achieved 
in deep optical lattices and with atoms having large vec- 
tor shifts. It is combined with the probability to have 
"inactive" atoms unaffected (P) pa 0.87. We give an er- 
ror summary for the no-blockade and blockaded cases in 
Table I. 

It was shown in |26[ that fault-tolerant MBQC can 
be realized with a 3D cluster state using topologically 
protected gates. In particular, an error threshold for an 
entangling two-qubit gate during the preparation of the 
cluster state has been calculated as e ps 7.5 • 10 -3 . Re- 
cently the threshold was increased to 1.1 — 1.4% j29j . Un- 
fortunately, we cannot directly compare the gate errors 
in our scheme with this threshold, since in [26|, [2{| the 
error was assumed to be partially depolarizing with the 
phase gate acting as Upg = (l — e)[I a Ib] + (e/15)([I a Xb] + 
... + [Z a Zb]). In our case the errors are due to leakage out 
of the computational subspace. We note, however, that 
the errors of the phase gate without and with blockade 
are higher than the threshold, but not too far from it. 

Our analysis shows that the no-blockade and block- 
aded case have similar phase gate errors, but it requires 
a very deep lattice and atoms with large hyperfine num- 
bers in the latter case. Additionally, the probabilty to 
have "inactive" atoms unaffected is higher in the no- 
blockade compared to the blockaded case. On the other 
hand, the optimal error combination in the no-blockade 
case requires a four-photon excitation, while in the block- 
aded case one-photon excitation gives the optimal er- 
ror combination. Generation of a ID cluster state in 
the no-blockade case requires almost an order of magni- 
tude less time (80 /is) compared to the blockaded case 
(700 fjs) due to the lattice change involved in the lat- 
ter case. The times required to produce a 2D cluster 
state are, however, comparable: 1.62 ms without and 3.4 
ms with dipole blockade. From the overall experimental 
complexity (deep double-well lattice, lattice manipula- 



tion required) the blockaded scheme seems more difficult 
to realize. 

The same scheme can be applied to generate the clus- 
ter state with polar molecules. For example, molecular 
states with small dipole moments can be used to encode 
a qubit while for the phase gate molecules can be ex- 
cited to a state with a large dipole moment, such that 
in this state molecules can interact via dipole-dipole in- 
teraction. A good candidate is 13 CO [3(|, where the 
ground electronic state has a rather small permanent 
dipole moment pa 0.1 D and long-lived nuclear spin sub- 
levels can be used to encode a qubit. It also has a 
metastable a 3 Hq state with a significant permanent 
dipole moment pa 1.4 D, in which molecules can inter- 
act via dipole-dipole interaction. Polar molecules, that 
currently or in the near future can be cooled to ultra- 
cold temperatures and placed in an optical lattice, are 
limited to alkali di-atoms, having permanent dipole mo- 
ments of (j, < 6 Debye (1 Debye=10 -18 esu cm). It lim- 
its the dipole-dipole interaction strength [i 2 / R 3 ~ 100 
kHz between nearest neighbors with R ~ 500 nm. The 
rather small interaction strength makes the no-blockade 
phase gate preferable for polar molecules. As an exam- 
ple, we can estimate the gate error for 13 CO. Assum- 
ing two-photon excitation to the large dipole moment 
state, the error due to the variation of the Rabi fre- 
quency e n var = 2tt 2 ^((50/(0)) 2 ^ ps 1.1 • 10~ 2 , with the 

lattice depth Vq = 100 E p. The error due to the finite 
ratio of the interaction strength to the Rabi frequency 
eimpoxc = Vgt/8(Jl a /A) a « 8 • 10" 5 assuming p = 1.4 
Debye, R = 500 nm, resulting in Vj nt ~ 2.5 kHz and 
fi 2 /A ~ 100 kHz. The large dipole moment molecular 
state is metastable with the lifetime ~ 500 ms, giving 
7 = 2 s" 1 and the decay- induced errors eR y db decay l — 
2irj/V int ps 8 • 10- 4 , e Rydbdccay2 = ^ 7 /(^ 2 /A) ps 10" 5 , 
resulting in the total gate error e ps 1.14- 10~ 2 . As we dis- 
cussed above for atoms, the gate pulses affect molecules 
in "inactive" sites if two-photon excitation is used. The 
probability to have the "inactive" molecules unaffected 
can be increased using multi-photon, e.g. four-photon 
excitation. In this case the error due to the Rabi fre- 
quency variation is four times larger, and the total gate 
error becomes e ps 4.3 • 10~ 2 , while the probability for 
"inactive" atoms to stay in the initial qubit states be- 
comes (P) ps 0.994. This analysis shows that gate errors 
for molecules and for atoms are comparable. 

In conclusion, we propose and analyze generation of 
a cluster state for measurement-based quantum comput- 
ing with neutral atoms and polar molecules in an opti- 
cal lattice using van der Waals and dipole-dipole inter- 
actions. We consider two schemes for implementation 
of a phase gate between pairs of nearest neighbors re- 
quired to generate the cluster state: without and with 
individual addressing within a pair. We show that in the 
former case the gate error e sa 1.25 • 10~ 2 is feasible with 
two-photon excitation to Rydberg states, provided the 
excitation-deexcitation to interacting states is adiabatic 
with respect to motional frequency of the lattice. Two- 
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TABLE I: Main errors 


in the no-blockade and blockaded cases. 


no — blockade 


two — photon excitation 


four — photon excitation 


Efiva, = (27T 2 ) ((6n/{ty)' Z ) « 1.05 ■ 10" 2 

timpcxc = 1/8 (Vi nt /n 2 /A) 2 « 1.25 ■ HT 3 

ERydb dccayl = 27T7/Vi n t ~ 6.7 ■ 10" 4 
ERydb dccay2 = ITy / (0? / A) « 3 ■ 10 _S 
Enon adiab = (tT 2 /4) ( ffl /i?) 2 « 2.5 ■ 10~ 2 

e « 1.25 • 10 -2 f 


En var = (8tt 2 ) ((SQ/{Q))' z ) « 4.2 • 10" 2 

eimp cxc = 1/8 (Vint/fi 4 /A 3 ) 2 « 1.25 ■ 10~ 3 
ERydb dccayl = 27T7/Vint « 6.7 ■ 10~ 4 
ERydb dccay2 = 7T 7 /(n 4 /A 3 ) « 3 ■ KT 6 
Enon adiab = (tT 2 /4) (a/i?) 2 « 2.5 ■ lO" 2 

E « 4.4 • 10" 2 f 


Einact cxc — 1 — (-P) ~ 0.25 


Einact exc = 1 - (P) « 6.5 ■ 10" 3 


with blockade 


two — photon excitation 


one — photon excitation 


envar = (2tt 2 ) ((6Q/(Q))' z ) w 0.15 


eq var = (tt 2 /2) ((^/(fi)) 14 ) « 1.13 • 10~ 2 
E imp week = (« 2 /A) 2 / (2|A vcc| 2 ) » 5 • 10~ 3 

ERydb decay = 7^/(40) » 1.75 ■ 10" 2 

emf fluct ~ (SljTpg) 2 ~ 6 • 10" 4 
e « 3.44 • 10~ 2 f 


E imp block = (fi 2 /A) / (2|A vcc| 2 ) « 2 • 10~ 2 

ERydb decay = 77T7/(4fi 2 /A) » 4.38 ■ KT 2 

£mf fluct ~ (SljTpg) 2 ~ 6 • 10~ 4 
e « 0.21+ 


Einact exc = 1 — (-P) ~ 2 • 10 3 


Einact cxc — 1 — (P) ~ 0.13 



not including einact exc and e non adiab 



photon excitation, however, leads to high error due to 
excitation of atoms in "inactive" wells at the minima of 
the standing wave. The optimal combination of the gate 
error and the "inactive" atoms excitation error can be 
realized using four-photon excitation to Rydberg states, 
giving the gate error e w 4.4 ■ 10 -2 and the probabil- 
ity to have "inactive" atoms unaffected (P) « 0.994. In 
the second scheme individual addressing within a pair al- 
lows to implement the phase gate using Rydberg (dipole) 
blockade. Addressing of atoms within a pair can be re- 
alized in a polarization-gradient double-well lattice. The 
gate error in this case is e ps 0.21, but can be reduced 
to e ~ 3.44 • 10~ 2 using one-photon excitation to Ryd- 
berg states, atoms with large hyperfine quantum num- 
bers, such as cesium and holmium, and a very deep op- 
tical lattice. With one-photon excitation the probabil- 
ity of "inactive" atoms to stay in their initial state is 
(P) « 0.87. 

We also analyze the lattice manipulation required in 
the latter case to realize the phase gate between all 
nearest neighbors and show under which conditions the 
manipulation is adiabatic. The total time required to 
produce a ID cluster state T\u ~ 80 fis without and 
Tid ~ 700 /is with individual addressing in a pair, where 
parameters of 87 Rb were used for the estimate. 

Finally, we show how the scheme for ID cluster state 
generation can be extended to realize a universal 2D clus- 
ter state and estimate the total time required T2D ~ 1-62 
ms without and T2D ~ 3.4 ms with individual addressing 
in a pair. 
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VI. APPENDIX A. CALCULATION OF THE 
AVERAGED FIDELITY OF THE PHASE GATE 
WITHOUT INDIVIDUAL ADDRESSING 

We calculate the fidelity of the phase gate Upa aver- 
aged over all initial two-qubit states 1 00} , |01), 1 10} and 
111) as follows 



F 



1 



|(00|*7 PG |00)| 2 - 
-(10|1> PG |10)| 2 



- (01|C/p G |01)| 2 - 

- (ll|C/p G |ll)| 2 



where we project the final state after the phase gate 
^p G |eiE2) (ei,2 = 0,1) on the state, expected after 
the ideal phase gate 1 11} - 111), |01> -> - 101), 
1 10) -> -|10), and |00) |00). We assume that er- 
rors due to different mechanisms are small and can be 
analyzed separately and added together. We limit our 
analysis to intrinsic gate errors assuming that technical 
errors can be in principle eliminated. 



A. Error due to a finite ratio of the interaction 
strength to the two-photon Rabi frequency 

Fint/(fi 2 /A) 

In this subsection we calculate the errors due to im- 
perfect excitation from the 1 1 1 } to the \rr) state. In Sec- 
tion IIA we derived a system of equations for the am- 
plitudes a%x, a+ and a rr of the two-qubit states |11), 
|+) = (|lr) + |rl))/V^, and | rr) in the case when the ini- 
tial state is 1 1 1} (see Eq.lJTJ). The corresponding Hamil- 
tonian is 



/ 2f2 i 



H/h = 



2H| 
~~A~ 

A 





2 A 2 

A 

\/2niQ2 

A 







A 



(11) 



2Q- L 
A 
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Assuming fti = O2 = ft for the simplicity of the anal- 
ysis and VS n t <C ft 2 /A we can calculate the eigenvalues 
of the Hamiltonian: 



1 



= 2Q Vint 

1 A 2 32 (ft 2 / A) 2 ' 



A 2 : 

4ft 2 



1, 



5 VSt 
64 J] 2 /A : 



5 K„t 



1 

4 int 64ft 2 / A' 



the — 1 11) state. The wavefunction at the end of the gate 
is as follows 



l*> = 



Vint ■• y int / Vint \ 2 

ai|*i)e 2l ^+a2|^2)e 4 "^Ta KcFJZ) 

■• Vjnt | 57ri / Vjnt \ 2 

+a 3 |* 3 )e l4 ^^M + 64 15*/aJ + 



1 ^ 1 / Vint 

2 I 16 \n 2 /A / 



1 / V 



int 



4 I 8ft 2 /A 16 V^7 A , 



sv/a _|_ 



+ 



with the corresponding eigenfunctions 



V2 I 32 V^ 2 /A 



1 / x _ 3 /- Vfat 



V2 \ 32 V^ 2 /A 



2 V 16ft 2 /A 512 (ft 2 /A) 
1 ( x 1 17 V int 



^2 V 16ft 2 /A 512(ft 2 /A) 2 




_1 / , 5 , \\ 



2 

int 



2 V 16ft 2 /A 512 (ft 2 /A) 2 

,,,, 4fi4i^iii., 



+ 



2 V 16ft 2 /A 512 (ft 2 /A) 5 



1 Vfat 17 V 2 



V2 V 16ft 2 /A 512 (ft 2 /A) 2 
5 Vint 23 V 2 



2 V 16ft 2 /A 512 (ft 2 /A) 2 

up to the second order in VS nt /(f2 2 /A). 
The state 111) evolves as 



tf) = ax 1*0 e- lXlt + a 2 |* a ) e~' iA2t + a 3 |* 3 > e 



iA 2 t 



-iA 3 i 



where the coefficients are determined from the condi- 
tion = 0)) = 1 11) and are given by: 

1 /, 1 / Vint 

ai = — | 1 



^\ 32 V^ 2 /A 

1 (. 3 Vint 25 f V^t 

a 2 = -- 1 



2 I 16ft 2 /A 512 V^ 2 /A 
1 (. , 3 V„t 25 ( Vmt 

«3 = - I 1 



2 I 16ft 2 /A 512 V^ 2 /A 



1 ( 1 Vint 1 / Vint ' 

4 I 8ft 2 /A 16\n 2 /A / 




y/2 \ 16ft 2 /A 512 V^ 2 /A^ 



*i)e 4 " 2 /a + 



1*2) x 



xc 



Vint 5 



Tri / Vint V 

2» Kn^JK) 




5 V n 



23 / Vi 



y/2 \ 16ft 2 /A 512 V^ 2 /A 



1*3) x 



i* Vmt 1 B»< ( Vint V 
X g J 8 C!2/A + 128 Vri2/A/ 



Expanding the exponents in Vint /(ft 2 /A) and keeping 
terms up to the second order we get 

[ 1 ; \A 128J \Q 2 /AJ 8 ft 2 /A 



with the corrc- 



As a result, |-(11|*)| 2 = 1 - J]*)* ' 

V- 2 

sponding error 2 (o^'/a) 2 • This error is present only for 
the initial state |11), since for the initial states 1 01) , 1 10) 
and 1 00) as is seen from Eqs.(3),(5) the state \rr) is not 
populated and no interaction is involved. The averaged 



over all initial states error is therefore 



This analysis also allows to calculate the error 
due to the unwanted interaction of atoms belong- 
ing to different excited pairs. If we assume that 
VntTint = 7r(l + <5Vnt/Vmt), where £Vint is due 
to the interaction between closest pairs, the cor- 
responding error averaged over all initial states is 
(3^ 2 /16) (1/8 + 19^ 2 /256) (<5V„t/Vi„t)(Vi„ t /(ft 2 /A)) 3 . 



The phase gate is realized in three steps: first, a it pulse 
is applied to both atoms during the time T such that 
ft 2 T/A = 7r/2. Next, atoms interact in the \rr) state 
and accumulate the phase Vint Tint = it. Finally, a second 
7r pulse of the same duration T de-excites atoms back to 



B. Error due to decay of Rydberg states 

In this subsection we calculate the error due to decay of 
Rydberg states. We can calculate this error using Eq.([5]) 
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for the |11) initial state assuming that the |r) state decays 
with the rate 7. Then, the |11) state at the end of the 
gate becomes 



I*) 



11 e 



-77r/(n 2 /A)-2 7 7 I -/Vi„ t 



111) 1- 



2771- \ 

n 2 /A ~ v~ t ) ' 



77T 



where we expanded the exponent using the smallness of 



2-fTT 



7 /(r! 2 /A), 7/^. As a result, |-(ll|#)r = 1 
fg, giving the error fife + fe 

If the initial state is 1 01} the system evolves ac- 
cording to Eq.Q. At the end of the gate I*) = 
- |01) e -T 7r /( 2 ° 2 /A)-7^/Vi„ t _ The 0V e r lap with the desired 
-|01) state is then |-(01|*)| 2 = e -7V(o 2 /A)-2 7 ,r/V, nt ) 
resulting in the error 771"/ (fi 2 /A) + 2771-/ Vint- Similar re- 
sult is obtained for the 1 10) state. There is no error due 
to the Rydberg state decay for the |00) initial state. The 
averaged over all initial states error is according to Eq. ^ 
given by 7 7r/(n 2 /A) + 2771-/^. 



Error due to non-adiabatic excitation to the \rr) 
state 



When two atoms in the |11) internal state and in the 
ground motional states of the lattice \g1g2) are trans- 
ferred to a doubly excited Rydberg state \rr), the state 
\rr) <X> \g1g2) is no longer the eigenstate of the total 
Hamiltonian, including vdW or dipole-dipole interaction. 
The interaction admixes higher-energy motional states 
to the original ground state as well as other Rydberg 
states to a much smaller extent. As a result, new eigen- 
states are superpositions of several motional states. Let 
us call the new ground motional state of two atoms in 
the presence of vdW or dipole-dipole interactions \g), 
while \&k) are the higher-energy motional states. If 
the interaction is not very strong only closest in energy 

states arc admixed, 



\~9)*\9i92) + a(\ 9l eP) + \eWg 2 )) 



in the first order perturbation theory, where 



is 



the j excited 

a = (5152 I Vint 



Si 4^ 



state of the 

■ E, 



i th atom, and 



J 9i 



(5lff2| V^nt 



(1' 



-.(1) 



92 ) 



E (i) ) (assuming 

e 2 

If the atoms are 



motional 

= {gm\VirA 

excited to \rr) adiabatically with respect to the energy 
separation between \g) and higher-energy states (which is 
of the order of the motional energy splitting) , the state 
l<7iS2) gradually evolves into \g) and there is no uncer- 
tainty in the interaction strength. If, however, the exci- 
tation is non-adiabatic, the state \g\g 2 ) does not change 
during the excitation, i.e. it is now a superposition of 
the new eigenstates \gig 2 ) ~ \g) + ft |ei), where (3 ps —a. 
As a result, after the first tt pulse 



111) 



\9i92) 



\rr) ® \ gi g 2 ) = \rr) <8 [\g) + j3 |ei)). 



During the interaction time the state |ei) acquires a 
phase factor exp(— iAETint/fi), where AE is the energy 
difference between the |ei) and \g) states. After the in- 
teraction time the system is in the state 

\rr) ® (\~g) + Pe-™^'^ 71 ^)). 

The second n pulse brings the system back to the |11) 
state without changing the motional state. As a result, 
after the second 7r pulse the state looks as follows 

= - |11) ® (\g l92 ) +/3(e— AE / y ^ n - 1) (ex)). 

The projection to the ideal — |11) ® \g192) is then given 
by l + \p\ 2 ( e - i7rAE ^ h -l) ^l + \a\ 2 (-iTrAE/V int h- 
(TT 2 /2)(AE/V int h) 2 ). Since AE « hu, and 

(si32|Vi„t gie[ 1] ^ ~ (a/R) {g x g 2 1 V^t \g192) = (a/R)V int , 

the perturbation theory gives \a\ ~ (a/R)(V{ nt /uj). As 
a result, the projection to the ideal final state is 1 + 
(a/R) 2 (V int /oj) 2 (-inoj/V int - TT 2 ui 2 /2V 2 t ), which gives 
the error ir 2 (a/R) 2 — n 2 (a / R) 4 (Vi nt / lu) 2 . Since this error 
is present only for the |11) initial state the averaged error 
is (n 2 /4)(a/R) 2 - (n 2 /4)(a/ R) 4 (V int /uj) 2 . 



D. Error due to the variation of the Rabi frequency 

In this subsection we calculate the error caused by the 
variation of the Rabi frequency experienced by each atom 
due to the finite width of the ground motional state. If 
the initial state is |11), the effect of the Rabi frequency 
spread can be calculated from Eq.@, where the Rabi 
frequency ft = (ft) +Sfl, and we again assume Oi = SI2 — 
fl. Assuming also that for the two tt pulses (ft) 2 T/A = 
7r/2, the state at the end of the gate is given by 



j \J>^j = e 4,ni5n/{n)+2Tri(Si1/{n)) 2 

|ll)-V2|+) + |rr) 



111) 



+- 



111) 



V2\ 



-4-nisn/ (n) -2iri(sn/ (n)) 



a 4mi&Q,/ <n) +2m( Sfl/ (SI) ) 2 



+ 



_l_ ^ _l_ }_ e -2iTisn/{n)~TTi(sn/{n)) 2 _^ 



-2TriSn/(n)+iri(5£l/(n)) 2 . 
~ o I X 



|ll)-|rr) |H)- V2\+) + \rr) 



xe 



-2m5n/(n)—. 7ri((50/<n)) 2 



\n) + V2\+) 



= 2iri(5Sl/(n)+7ri(5n/(n)) 2 



Expanding the exponents up to the second order in 
8Q/(Q) we have |*> = - 111) e 4«<5n/<n)+2«(5si/<si» 2 ) 
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which gives the projection to the desired — |11) state 

Kii|*)| 2 ~o((«V(^)) 3 ). 

If the initial state is 1 01) the wavefunction is given by 
Eg.® 



|01) + \0r) ^t/A 



Assuming again O = (f2) + Sft and (f2) 2 T/A = tt after 
the two tt pulses, the wavefunction becomes 

|^\ _ _ e 4,niSn/{n)+27ri(5n/{n)) 2 -i7rA hf /A x 

x (|oi) cos (2nsn/(n) + Tr(sn/(ri)) 2 ) + 

+i |0r) sin (2tt<50/(Q) + 7r(<5f7/(f7>) 2 )) . 
The projection to the ideal — |01) state is then 
- (01|<£) = e-iT&hf / A+Airtsn/ (n)+2iri(sn/ (n)) 2 x 



x 1 - 2tH 



giving the error 4tt 2 (6il / (il)) 2 . A similar result is ob- 
tained for the 1 10) state. 

Finally, if the initial state is 1 00) , it evolves according 
to Eq.©. After the phase gate the state turns into | \1/) = 
1 00) e ^isn/(n)+t7:(sn/(n))\ The overlap with the desired 

|00) state is |(00|*)| 2 = 1. 
The averaged over all initial states according to Eq.© 

error is then 2tt 2 ^(<5f7 2 /(^)) 2 ), where the averaging m 

the error expression is over the ground motional state 
wavefunction. 



E. Excitation of atoms in minima of standing wave 
excitation pulse 




The overlap with the initial |11) state is then 

2Trin 2 /n 2 , , . . 

(11|*) = - [cos 2 [rrh 2 /^ -l + 2cos [TTtt 2 /VL- 

The initial 1 1 ) state (the same for 1 10) state), accord- 
ing to Eq.@ after the gate becomes 



|*) = e 2 ™?£ M01)cos ( tt^ j +i |0r)sin Itt^ 



which results in the overlap with the initial 1 01) state 
(01|tf) = e 2^ 2 /n 2 cos f^/Sl 2 ) _ 

The initial state |00) becomes |W) = |00) e 2l7r ° 2 /° 2 af- 
ter the gate, according to Eq.©. 

The Rabi frequency at the minima of the standing wave 
Cl = flo cos(37r/8), at the maxima SI = Jlo cos(77r/8), 
as a result, the ratio f2 2 /£l 2 0.17. The probabil- 
ity to find the pair of atoms in the initial state after 
the gate, averaged over all initial states, is then (P) = 

(i/4) (|(ii|*)| 2 + Koii*)! 2 + \(iom 2 + l(oo|*)| 2 ) « 

0.75. 



VII. APPENDIX B. CALCULATION OF THE 
AVERAGED FIDELITY OF THE PHASE GATE 
WITH INDIVIDUAL ADDRESSING 
(BLOCKADED CASE) 

Next, we calculate the errors of the phase gate with 
individual addressing in a pair of atoms. We start by 
calculating the error due to the finite ratio of the two- 
photon Rabi frequency to the vector shift A vec between 
atoms in the left and right wells and between the Rabi 
frequency and the interaction strength Vj n t, resulting in 
imperfect individual addressing. 



We also need to estimate the probability that the 
atoms in "inactive" lattice sites, i.e. at the minima of 
the standing wave excitation pulse, are not affected. Let 
us denote Cl 2 /A the two-photon Rabi frequency at these 
sites. Now, as can be seen from Eq.© the |11) initial 
state after the gate becomes 



I*) 



>1) 




+ 4 cos ^ 




( C ° S (^i?) 


4 







+ 



\rr) 



n 2 
'n 2 



iy/2\ 



sin 2tt 



Q 2 

v: 2 



o 2 
~n 2 



A. Error due to a finite ratio of the Rabi frequency 
to the vector shift (S1 2 /A)/A VGC and interaction 
strength (f2 2 /A)/Vi nt 

1. 1 11) initial state 

First we consider the |11) initial state. The phase gate 
is implemented in the following way: a tt pulse resonant 
to the |1) — \r) transition of the control, e.g. the left 
atom is applied. The states |11), |lr), |rl) and \rr) (other 
states are far detuned and have much smaller amplitudes) 
evolve according to the Hamiltonian: 
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( -(- 



\A+A V 



H/h 



£li£l 2 
' A+A v , 

Oi0 2 



' A+A v . 

n 2 



A+A v 



A+A v 



' A+A v . 



^ A+A vec 

OiQ 2 
A+A V1 



A 



Mr 





' A+A vc , 

' A+A vc , 
Q 2 
A+A v 



A+A v 



/ 

(12) 



We assume fii = fi 2 = ^, Mnt > A voc > ft 2 /A, 
and A ;g> A vcc , l^nt- In the analysis we will, there- 
fore, keep terms up to the first order in (i7 2 /A)/V;nt 
and up to the second order in (f2 2 /A)/A vec . Neglecting 



ai r ~ (^ 2 /AA vcc )aii and a rr ~ (i7 2 /AVint)a r i as small, 
we have a system of equations for an and a r i governed 
by the Hamiltonian: 



H/h 



_ ( 2f2 2 

4(1+ 



(n 2 /A) 2 

A vcc -20 a /A 
(Q 2 /A) 2 

2Q 2 /A) 



Vint (A 



40 



(tf/A) 2 

2H 2 /a; 



2tr 

A 



(^ 2 /A) 

Vi„t 



(13) 



The corresponding eigenstates and eigenfunctions are 
(n 2 /A) 2 ( 20 2 /A N 



After the first 7r pulse the wavefunction, therefore, turns 
into 



2n 2 

' A 



2A V 



1 



(» 2 /A) 2 ± O 2 
2A vnr A 



A 

(^ 2 /A) 
8A 2 CC 



vec 
2 / A \2 



and 
1*2) = 



1 
1 

'71 



1 - 



ft 2 /A n 2 /A i7(n 2 /Af 



1 2 

vcc 



4A VCC 4V int 32A! 
ft 2 /A ft 2 /A 15(fi 2 /A) 2 



V2 



1 



4A VCC 4V int 32A 2 C , 
ft 2 /A f7 2 /A 15(ft 2 /A) 



1 



4A VCC 4M„t 32A 2 CC 
ft 2 /A ft 2 /A 17(fi 2 /A) 2 



(14) 



111) 



Irl) 



111) 



The 



y/2 4A VCC ' 4Vi n t 
initial state 111) evolves as 



32A 2 ec 



I r 1^15) 



I*) = e 



Avoc 



2ii-i 2 t 1 i(iV/A)^t / , . 2S2f/A \ i{SVf/A)^t 



x ai |*i) e 



a 2 1^2) e 



I*) 



i7rf! z /A 
4A VCC 



ifi 2 /A 7rfi 2 /A ^ 2 /A i(fi 2 /A) 2 



2A V 



16A V 



A 2 

int 



Hi) 



+i 1 - - 



4 vec ^ ^int 

Next, a 27r pulse resonant to the target, i.e. the right 
atom is applied. The system of equations for the ampli- 
tudes an, ai r , a r i and CLfy IS the same as Eqs. (|12[ ) with 
A V cc — > — A V ec- The system of equations for an and ai r 
is also the same as Eqs.(fT3"|) with A vcc — > — A vcc . As a 
result, during the 2tt pulse the state evolves as 



2in 2 t 



x ai |*i) e V VC V^ 



+a 2 |* 2 )e 



— >• 



where 



01 



«2 



where 1^1,2) are given by Eqs. lfTS"!) with A v 
and |rl) — > |lr), ai,2 are given by Eqs.([T6|). After the 2ir 
pulse the states |11) and |rl) change as 



1 / n 2 /A fi 2 /A_17(OVA) 2 

V ~ 4 Avcc + 4Vi nt ~ 32A2 ec 

' n 2 /A n 2 /A 15(Q 2 /A) 2 \ 

V2 V 4A VCC 4Vl„ t 32A 2 CC J ' 1 j 



111) 



A 



1 



111) 



i7T /ft 2 /A 



8 V A v 



lr) 



16 



- A voc - Q 2 /A /, . 2Q 2 /a\ . Q 2 /A 

Irl) -> e 1 c V c / int | r i\ 

As a result, after the 27r pulse the state becomes 



of the states 1 01} and \0r) (other amplitudes are much 
smaller) evolve according to the Hamiltonian 



2Qi 



|tt> 



ix!2 2 /A I 1 2fi 2 /A \ iTtCl 2 

4A voc y-+ a voc J+— 



n 

' A+A h/ 



-Ah/ V 



O1O2 
A+A h/ 



7T 

16 



fi 2 /A zfi 2 /A + ^(fi 2 /A) 2 



A, 



2V in 



A 2 

^vcc 



' vec A+A h/ +A v 

The initial state 1 01) changes as 



A+A v 



(17) 



xe 



2 A V cc 



(, 2Q 2 /A\ ijrfi 2 
1 s^- J+— , 



|11). 



1 



7T 

512 



(ft 2 /A) 2 



xe 



n 2 /a a voc 



rf2 2 /A A 2S"i 2 /A \ _ ii 
A v cc I A vcc J 



CI 2 /A 



•1) 



where we neglected the small ^ A i A ^ | lr) term. 

vcc 

Finally, after the second n pulse resonant to the left 
atom, the states 1 1 1} and |rl) evolve as 



111) 
|rl) 



ai|*i)e iAlt + a 2 |* 2 ;e 
a 2 |*i)e- lAlt -oi |tf 2 )e _iA2 '. 



-i\ 2 t 



At the end of the phase gate the initial 1 1 1} state 
evolves into 



|*) -> -e 



2 /A 2S^ 2 /a \ 

/CC I ^VCC J 



rvii 2 /A 



Q 2 /A Q 2 /A 
1 — in— in- 



_ /33Z7T (» 2 /A) 

Ue + 2 J A 2 ec 

'n 2 /A n 2 /A 



VS„ t 
2 



111) 



A v , 

7T 2 

32 



37TZ 



(n 2 /Af 

A 2 



|rl> 



The projection on the desired — |11) state is then 



- (11|*) = e 



TV i'l 2 / A 



ft 2 /A 2 /A 
1 — in— in— — — 

^vec Mnt 

^33™ i n 2 \ (ft 2 / A) 



V 16 



A 2 



The error can be found from ( — (111*)! 
0{(n 2 /A) 2 /A vcc V int ). 



1 + 



2. 1 01) initial state 

Now we analyze the evolution of the 1 01) state. During 
the first n pulse, resonant to the left atom, the amplitudes 



I*) 



x l*+>e 



hi x 

Avoc' 



t |.-ffl||..). 



2 A 2 

^ ^vec 

where the eigenstates are given by 

n 2 /A 



oc * ( 1 „ (o 2 /a) 2 "\ ' 



l*+> 



01 



l-^'llOr) 



A 2 



VCC 

2, 



(18) 



After the first n pulse the state turns into 

|*) ->• ie l2<A +^/' 



1+2 ; 7 + 



+(e --A vec / ( nVA) _ 1 _ ^) ("V A)2 > ) |01) + 



+- 



ft 2 /A / 7rft 2 /A 
-r^— 1 + 2 , 



2A V 



VCC 
rrA v 



I - ^ /A ] e-TO ) |0r) 



During the 27r pulse resonant to the right atom, the 
1 01) and | Or) states evolve as 



|01) 



|0r) 



e' 4 + 4 '/ +1 A - ((|01) - |0r))e-*° 2t / A + 
+ (|01) + |0r))e lo2 */ A ) , 

e A+A */ + A - ((|01) + |0r))e jO t/A 
— (jOl) - |0r))e-^ 2 */ A ) . 



As a result, after the 2n pulse 

|01) -> |01)e i7rA/(A+A ^ ) , 
\0r) -> |0r) e" A /( A + A "/), 

During the second n pulse the states 1 01} and |0r) 
change in the following way 

-Avcc I ■ iA _ - irfl 2 /A 
"> 2 / A A + A i, , 4 ^t; 



|01) -> *a+|tf+)e 'i^^+^z '^- + 

n-n 2 /A 



| Or) 



+ia- \-qj_) e iA+A ^ +l ~— , 
io_|*+)e 42 ^7A+ 4 A+A ft/ 

ttA i - 7rrZ 2 /A 

-o+ |*_) e J ^T7+ l ^c- ; 



- 7Tf2 Z /A 

4 — A 
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where \^±) are given by Eqs. (|T8| and a+ 

-(ft 2 /A)/A vcc , a_ = 1 



(SI /A) 
2A?„ „ 



After the second 7r 



pulse the initial 1 01) state becomes 

_ft 2 /A 



-2mA hf /A 



^vec 

e __ 2 _^« )l01) 



A 2 

vec 



+17T 



» 2 /A 

"A 

(fi 2 /A) 



1 - e + 



vec 
2 / A^2 



A 2 

vec 



i + P r"^M 



I Or) 



The projection on the — 1 01) state is then 

- (01|1<) = e - 2 «A h/ /A + 
+ e - 2 



D 2 /A 



A 

7T_ 2 ^ (0 2 /A) 

"2 



vec 

2\ ( '02/A\2 



A 2 



Assuming 7rA V0C /(rj 2 /A) = 2?m we have |-(01|*)| 2 = 
1 — 2 ^° A 2 A ^ , and the resulting error is 2 ^° A { A ^ . 



3. 1 10) initial state 



much smaller) change according to the Hamiltonian 



H/h 



A+A 



niii 2 



A 

V ^ A+A h J 



(19) 



During the first 7r pulse the states |01) and |0r) evolve as 



|10) 



|r0) 



-(|10> + |r0»e' 



.- n 2 t i -■ n' 2 t 1 

e ls ^ + ^ 5 ((|10) + |r0)y 



-(1 10) - |rQ))e- 



After the first ir pulse the states 1 10) , |r0) become 



|1Q) 



-ie 1 2A | r Q) , 



\r0) -> -ie~ l Ts-|10) 



During the first 7r pulse the amplitudes aio and a r o 
(other states are far detuned and their amplitudes are 



During the 2ir pulse resonant to the right atom, the 
amplitudes aio and a r Q are governed by the Hamiltonian 



J 



H/h = 



A-A vcc ^ A+A hf 
A-A voc 



OiS7 2 
A-A vcc 

A+A hf ~A v 



A-A v 



(20) 



As a result, the state |r0) evolves as 

Q 2 /A 



|r0) 



A v 



where 



|*+) 



n 2 /A 



A 2 



1 10) 



A \~ A 2 

i-ivec V '-Vec 



A vec 

(n 2 /A) 2 



|r0) 



After the second ir pulse the initial state 1 10) evolves 
into 



I*) 



-2-rri- 



1 + ITT 



1 - e "Ta l + 



» 2 /A 
"A 

(» 2 /A) 
A 2 

vec 



vec 
2 /A^2 



1 10) 



+ e " 2 / A I 1 + ITT 



^vec \ ^vec 

^ 2 /A / n 2 /A\\ . . 

^vec \ *-*vee / / 

The overlap with the desired — 1 10) state is then 



(io|*> = 



\ A 2 

\ vec 

/ n 2 /A ( ir 2 \ (n 2 /A) 2 ' 
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If ttA V0C /(07A) = 2nn 



- (10|#) = e 



2 A 2 



resulting in |-(10|*)| 2 = 1 + 

Finally, the 1 00) state evolves into |00) -> 
| 00 ^ e 2mA/(A+A ft/ ) after the gate and |(00|*)| 2 = 1. 

The error averaged over all initial states is then e = 

(n 2 /A) 2 



B. Error due to the decay of Rydberg states 

We again start the analysis from the initial state |11). 
After the first tt pulse resonant to the left atom it be- 
comes 

|*) -> _ ze -7-/(20 2 /A) | rl) _ 

At the end of the 2tt pulse resonant to the right atom 
the state |rl) is given by 



Irl) 



^A vec /(n 2 /A)- 7 7r/(n 2 /A) 



|rl). 



Finally, after the second tt pulse resonant to the left atom, 
|rl) _ 2e -7/(2^ 2 /A) | U) _ 

As a result, at the end of the gate 

_ e ^A voc /(0 2 /A)-27r 7 /(0 2 /A) _ 

From the overlap | — (11|\I/)| 2 = e - 47r f/(^ 2 /A) the error is 

47T 7 /(fi 2 /A). 

The initial state |^) = 1 01) after the first tt pulse 
turns into |*) -> -|01). After the 27r pulse |*) -> 
e -inA hf /A | Q1 ^ Finally, after the second tt pulse |#) 

_ e -i7rA h/ /A | Q1 ^ giying |_(01|<If)| 2 = 1. 

The initial state = 1 10) after the first tt pulse 
evolves into |W) — > — 

W(2n 2 /A) | r0 ^ After the 

2tt pulse |r0) -> e -^A h/ /A+i7rA voc /(f2 2 /A)- 7 7r/(Q 2 /A) | r Q^ 

The state |r0) turns into |r0) -> -j e - i7rA "// A 1 10) af- 
ter the second tt pulse. As a result, at the of the 



gate |*) -> - e -2^A h/ /A+i7rA voc /(f2 2 /A)-37r 7 /(20 2 /A) hQ\ 

This gives |-(10|*)| 2 = e -3T7/(" 2 /A) reS ulting in the er- 
ror 3nj/(n 2 /A). Finally, the initial state |00) -> |00) 
after the gate and there is no error due to the Rydberg 
state decay. 

The error averaged over all initial qubit states is then 
77r 7 /(4fi 2 /A). 

C. Error due to the variation of the Rabi frequency 

1. 1 11) initial state 

First, we analyze the |11) initial state. During the first 
tt pulse, resonant to the left atom, the state evolves as 

|^ = e 2tn 2 t/A ^ jjjj - \ rl ) c -m 2 t/A + I 11 ) + \ rl ) c in 2 t/A^j 

Again, we assume that fl = (fl)+6Q. and (il) 2 T/A = tt/2 
during a tt pulse. As a result, after the first tt pulse 

1^ _?_ e 2Tn5n/(n)+i7r(m) 2 /(n) 2 x 



2ni-^- + iTT^X- 




2-7T 





2 (sn) 2 



|H) 
|rl) 



During the 2tt pulse resonant to the right atom, the states 
1 1 1) and |rl) evolve as 

|ii) -» e 2in2 */ A /^L!±L_L!a2, 

+ 2 

|rl) — > |rl) g iA voct-2Tri0 2 t/A^ 

so that at the end of the 27T pulse 



+ 



111)^-^ 



f p 47TiSn/(Q)+2TTi(SQ) 2 /(Q) 2 x 



2 _ 47r 2(M! 



I |lr> , 

The initial state |11) after the 2tt pulse becomes 



( - ( 2 -$y 



+ Tri 



(sny 



l*> 



67ri6fi/(fi)+3i7r(5fi) 2 /(fi) 2 



e4 7rA voc /(^ 2 /A)-4 7rtI 5JV(^>-27r 4 («2) 2 /(f2) 2 ^ - 8ni-^ - (4ni + 17tt 2 ) 



2 n(50) 2 



<J2) 2 
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After the second ir pulse resonant to the left atom the 
states |11), |rl) and |lr) become 



|H) 
\rl) 



2 ml 
1 m 2 



in) 

|rl) 



|lr) ->• |lr) e -^ A voc/(0 2 /A)-47r 4 50/{0)-27 rt ( I 5^) 2 /(0> 



If 7rA vcc /(r2 2 /A) = 27m the wavefunction after the 2w 
pulse is given by 



' -§ (a - 16tt^ - (8tu + 32tt 2 ) 



(<5q) 2 
W 



r 



I*) 



1 c 2m.5n/<n)+i7r(<5n) 2 /<^> 2 v 
'2 6 



47r "(Hp" 



-16- 2 ^ ) 




(21) 



The projection to the desired — |11) state is then 
- (11|*) = e 2 * iSn /( n )+™( sn ) 2 /( n ) 2 x 

X I 1 — 4-7T2- 



sn , n n 2 ,((50) 2 



giving |-(11|*)| 2 



2. |01) initial state 

Next, the initial |01) state during the first ir pulse 
evolves as 



|1>) 6 ^ 4 a -a- 



|01) 



After the ir pulse 



I*) 



101) 



During the 2tt pulse resonant to the right atom, the 
1 01) state changes as 



2 * & h f 1 



|01) — > e 2j — — 



((|01) - |0r))e- tO2 */ A + 
+ (|01) + |0r))e iO2t / A ) . 



After the 27r pulse the initial state becomes 



If) 



<n> 2 







27T7-S- + ?7T^^ 



( ioi> ^ 

l |0r) ) 



(22) 



During the second tt pulse resonant to the left atom, the 
states 1 01) and |0r) evolve as 



|01) -» e * * a a |01) 



|0r) 



iA vec t+2i 



|0r>. 



As a result, at the end of the gate the initial state 1 01) 
turns into 



1 - 2tt 



l*>- 

2 (<50) 2 







2 ™w + i7r w 



From 



-<oi|*)r 

4^ 2 (<50) 2 /(0) 2 . 



47T 



2 (sny 



( m ) 

{ |0r) J 



the error is 



3. 1 10} initial state 
The initial state 1 10) during the first ir pulse evolves as 

I*) = e 2^"^^ 1 ( ( | 10 ) _ | r0))e -^*/A + 

+ (|10) + |r0))e^ 2t / A ) . 
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After the tt pulse 



I*) 



2ni™+in^ 




o \ / |10> \ 



D. Excitation of atoms in minima of standing wave 
excitation pulse 

Finally, we analyze the undesirable excitation of atoms 
in "inactive" wells, situated at the minima of the standing 
wave excitation pulse. 



During the 2ir pulse resonant to the right atom the states 
1 10) and |r0) change as 



|10)^ e 2 ^-V^|io), 



iA„ ec +2i 



|r0) -> < 



i.e. after the 27T pulse 



^ |r0> 



l-ir>\ 4 ^Tm+ 2 ^ iL 77TTr 

1 10) -> e <!2> <"> 2 

-T^T+4^+2™^ 



|10) 



. Ah/ 



|r0> e ^ w 3 " |r0) . 

At the end of the 27r pulse the initial state becomes 



27rt w +»7r^ 



\ / |10) 



After the second 7r pulse the states 1 10) and |r0) turn 
into 



( |10> ^ 
I |r0) ) 



<n>-* 



'<n) 



.(so.) 1 



2-tt 2 & 2tt^ - /tt^t 



o _2 (<5^) 2 
/_7F "(OF 

'(12) ^ ra 7?2p^ 



|10) 
|r0) 



If 7rA V0C /(f7 2 /A) = 2im, the state |10) at the end of the 
gate becomes 



I*) -> -2 e ""'" r " 



1 9tt 2 ^ 



<5fi 



r 2(«2)f 

<n> 2 



|10) 
|r0) 



From |-(10|*)| 2 = 1 - 4tt 2 (Sty 2 / (fl) 2 the error is 
4^ 2 («}) 2 /(0) 2 . 

Finally, the 1 00) state at the end of the gate be- 
comes |00) -> e 8™«VW+4™(«2) 2 /<J2) 2 -4™A h/ /A | 0Q ^ giy _ 

ing|(00|vl/)| 2 = l. 

The error due to the variation of the Rabi fre- 
quency averaged over all initial two-qubit states is then 
27T 2 ^(<55!) 2 /(Sl) 2 ), where the averaging in the error ex- 
pression is over the ground motional atomic state. 



1. 1 11) initial state 



The initial state |11) after the first 7r pulse turns into 



+ 



|#) -> — — ((|11) - \rl))e-™ n l 2ir 

+ (|ll) + |rl))e" rA2 / 2f22 ) , 



where ft and Q are the Rabi frequencies at the 
lattice sites, corresponding to the minima and 
the maxima of the standing wave, and Cl/Cl = 

(l - v/(l + Vi/4T/ )/2) / (l + v/(l + V 1 /4V )/2) . Af- 
ter the 2n pulse resonant to the right atom, the states 
1 11) and |rl) evolve as 

|11) -> e 2 ^ 2 ' n2 x 
x (cos (irn 2 /^ |11) +i sin (irCl 2 /^ |rl>) , 
|rl) e ™ A vco/(n 2 /A)-27r J f2 2 /^ 2 | rl ^ _ 

During the second 7r pulse, resonant to the right atom, 
the states |11), |rl) and |lr) change as 

|H) 
|rl) 

cos (tt^ 2 /^ 2 ) ism^n 2 /^ \ / | U ) 
zsin(^ 2 /r! 2 ) cos (tt^ 2 /^ 2 ) J V l rl ) 

At the end of the gate the state |11) becomes 

I*) = e 4«^ 2 /o 2 cos 2 ^2 /2n 2^ cos ^ 2 /n 2 ) |ii) 

_ e - A vec/(n 2 /A) gin 2 ^q2 /2 ^ j n) + 
sin 2 (2tt0 2 /0 2 ) |rl) + 



+ 



i i7rA vcc /(J2 2 /A) . . 

+ sin [ntf/Sl 2 } |rl) + 



+ie 

x 



2 

27rif2 2 /fi 2 -i7rA vcc /(2fi 2 /A) 



cos (^7rSl 2 /2^ 2 ^) sin (7rtt 2 /0 2 ) \lr) . 
As a result, the overlap with the initial |11) state is 

(ii|*) = e ^ 2 /o 2 cos 2 ^n2 /2n 2^ C0S ^Ci 2 /n 2 ) 



_ e i7r A voc /(f2 2 /A) 



sin 2 (iTn 2 /2n 2 
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The probability that after the gate the atoms stay in |11) x (cos (tt0 2 /20 2 ) |10) + i sin ^7rf2 2 /2Sl 2 ^) |r0)) . 

18 During the 2n pulse, resonant to the right atom, the 

|(11|*)| 2 = cos 4 (.O 2 ^ 2 ) cos 2 (.0 2 /0 2 ) + sin 4 (.O 2 ^ 2 ^ |10> and |r0> chaDge &S 
-i sin 2 (ttOVO 2 ) cos (tt0 2 /0 2 ) cos (47r0 2 /0 2 ) , |10) -> < ' 



jtvQ 2 /n 2 +i7r(n 2 /n 2 )(A/(A+A h/ )) 



where we assumed 7rA vcc /(0 2 /A) = 27rn. 

The ratio O/O w 0.117 for Vi = Vb, giving |(11|*)| 2 
0.996. 



2. |01) initial state 
If the initial state is |01), it evolves into 

_ e i7r^ 2 /20 2 +i7rf2 2 /(20 2 )(A/A+A h/ ) |q-|\ 

after the first tt pulse resonant to the left atom. After the 
2n pulse resonant to the right atom the 1 01) state turns 
into 

| 01 ) _^ e nrn 2 /n 2 +tK(n 2 /n 2 )(A/(A+A hf )) x 
x (cos (tt0 2 /0 2 ) |01) +isin (7r0 2 /0 2 ) |0r)) . 

During the second tt pulse resonant to the left atom, the 
states |01) and |r0) change as 

|01) e ^^ 2 /2^ 2 +^(0 2 /2^ 2 )(A/(A+A h/ )) | Q1 ^ 
i7i-A vcc /(20 2 /A)+i7r0 2 /2^ 2 +i7r(0 2 /2^ 2 )(A/(A+A h/ )) 



|io), 

|r0) ->• e i7rA voc/(^ 2 /A)+i7rSl 2 /^ 2 +i 7r (0 2 /0 2 )(A/(A+A h/ )) | r Q^ _ 

During the second n pulse the states 1 10) and |r0) evolve 



as 



(|10) 
[\r0) 



e i7r0 2 /2f2 2 +i7r(n 2 /2f2 2 )(A/A+A h/ ) x 



cos (tt0 2 /20 2 ) isin(7r0 2 /20 2 ) \ / | 10 ) \ 
JS in(7r0 2 /20 2 ) cos (tt0 2 /20 2 ) ) V l r0 > J 

At the end of the gate the initial state |10) becomes 

= 27ri0 2 /0 2 +27ri(0 2 /0 2 )(A/(A+A h/ )) 



l*> = e 2 



x cos 2 (7r0 2 /20 2 ) 1 10) 



1 ,„ 



_ j_ i7rA voc /(a 2 /A)+27i-if2 2 /n 2 +2m(n 2 /n 2 )(A/(A+A ft/ )) 

4 



|0r) e 

At the end of the gate 

| \I>) = e""""" '"" 1 '"" M "' v " 1 "" J " x 
x cos (tt0 2 /0 2 ) |01) + 

+ie -iTrA vcc /(n 2 /A)+2iTTU 2 /n 2 +2iTT(U 2 /0 2 )(A/(A+A h/ )) x 

x sin (tt0 2 /0 2 ) |0r) . 



|0r). 



x sin 2 (tt0 2 /20 2 ) 1 10) + 

Vo 2 +27ri(f2 2 /f2 2 )(A/A+A h/ ) x 
sin (^0 2 /0 2 ) (l + e -A voc /(o 2 /A)^ | r0) ; 



which gives the overlap with the initial 1 10) state 



(101*) = e 2i ^ 2 /2^ 2 +2^i(0 2 /0 2 )(A/A+A h/ ) x 

(cos 2 (^0 2 /20 2 ) - e -A voc /( f 2 2 /A) sin 2 (^ 2 /20 2 )) , 



The overlap with the initial 1 01) state is then (01|\I/) = 

e ^A 2 /fi 2 +2 47 r(A 2 /^ 2 )(A/(A+A h/ )) CQS (^} 2 /0 2 ), resulting 

in the probability to find the atomic pair in the initial 

state after the gate |(01|*)| 2 = cos 2 (?r0 2 /0 2 ^ s=s 0.998. resulting in the probability to find the pair of atoms in 

the initial state |(10|*)| 2 = cos 2 (tt0 2 /0 2 ) , where we 
assumed that 7rA vcc /(0 2 /A) = 2nn. 

The state 1 00) after the gate bcomes I*) = 
The initial state 1 10) during the first tt pulse resonant |00) e 4 ^(" 2 /^ 2 )(A/A+A h/ ) ; giying |( 0|*)| 2 = 1. The 

probability to find the pair of atoms in the initial state 
after the gate, averaged over all four initial states, is 
(P) = 0.998. 



3. 1 10) initial state 



to the left atom evolves into 

|^ _ e ^0 2 /2^ 2 +i7r(f2 2 /2^ 2 )(A/(A+A h/ )) x 
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